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Abstract 

We analyze B — > KH C (where H c = rj c , J/ip, Xcj(j=o,i)) decays estimating non- 
factorizable contributions from the light-cone sum rules (LCSR). Nonfactorizable 
contributions are sizable for B — > KJ/t/j and particularly for B — > K\ c \ decay. For 
the B decays into a (pseudo) scalar charmonia we find small nonfactorizable con- 
tributions which cannot accommodate relatively large branching ratios obtained by 
measurements. This specially concerns the puzzling B — ► KxcO decay. 
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1 Introduction 

Last years there was a considerable progress in the measurements of B de- 
cays into diverse charmonium final states. These decays governed by a color- 
suppressed b — > c transition could provide a valuable information on the fac- 
torization properties of S-meson decays. 

The first observation of B~ — > K~Xco decay a year ago by Belle collaboration 
[1], Table 1, clearly shows the breakdown of the naive factorization assumption 
in the color-suppressed B decays into charmonium. Namely, this decay, and the 
corresponding B~ — > K~\c2 decay are forbidden in the factorization approach 
[2], due to the V-A structure of the weak vertex, i.e. 

(Xco\(cc) VtA \0) = (xc2\(cc) VtA \0) = 0. (1) 
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Table 1 

Experimental summary on branching ratios of B decays into charmonium [1]. The 
first result for a particular decay mode is Belle result, the second one is from BaBar 
Collaboration, and the third one (when exists) is from the CLEO experiment. 



decay mode 




B~ K- Vc 


12.5 ± l-4±i;§ ± 3.8 

15.0 ±1.9 ±1.5 ±4.6 
6.9+l;f ±0.8 ±2.0 


B~ -► K-J/i) 


10.1 ±0.3 ±0.8 
10.1 ± 0.3 ± 0.5 


B~ K- X co 


6.0±?;| ± 1.1 
2.7 ±0.7 


B- - K'xd 


6.1 ±0.6 ±0.6 
7.5 ±0.8 ±0.8 


B(B~ - ^- Xc0 ) 
6(5- -► K~J/ip) 


°- 60 to 18 ± °- 05 ± °- 08 



More surprisingly the branching ratio of B~ — > K~Xco is comparable with the 
branching ratios of the decays B~ — > K~J/tp and — > X _ x c i, which posses 
the nonvanishing factorizable amplitudes. Another nonfactorizable B decay 
into Xc2 charmonium was observed inclusively with a large branching fraction 

B(B - X XC 2) = (1.80±g;|| ± 0.26) • 10~ 3 , (2) 

which is of the order of the branching fraction of the factorizable B — > Xxci 
decay. 

However, even for the B decays into charmonium states which can be calcu- 
lated by the factorization assumption there is a problem of theoretically too 
low predictions which cannot accommodate the experimental data. Therefore, 
for all B decays into charmonium states one expects large nonfactorizable 
effects. 

In this letter we would like to approach the calculation of the nonfactoriz- 
able contributions to the B decay branching ratios into charmonium from the 
LCSR point of view. The method was developed for the calculation of the 
soft nonfactorizable corrections in B — > titt [3] and was later extended on the 
B — > K J I ip decay [4] . Here we will made a short presentation of the method, 
and the interested reader should have a look in the above references for the 
further details. 
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Table 2 

Charmonium states considered in the paper and their properties. 



H (J PC ) 

11 c\° J 


1 TT 


y- J \JHc\ rL cl 


mu ropviffii 

ln H c \y v v J L U J 


frr [MpVl 

J H c L ivi e v j 


/sTTfOeVl [71 


Vc(0~ + ) 




fr]c m ric 


3.0 


420 ± 50 [8] 


3.8 ±0.2 


J/V(l~) 




fj/i> m j/ip e n 


3.1 


405 ± 14 [4] 


3.8 ±0.2 


Xco(0 ++ ) 


cc 


fxcO m Xc() 


3.4 


360 [9] 


4.0 ±0.2 


Xcl(l ++ ) 






3.5 


335 [9] 


4.0 ±0.2 



The effective weak Hamiltonian relevant for our discussion on the b 
transition 



ccs 



Hw = ^[V^CMO, + C^Oi] - V tb V* £ Ci^)^} + h.c. (3) 
contain the current- current operators 

Oi = (cb) V - A (sc) V - A , 2 = (cc)v-A(sb) v _ A , (4) 

and the QCD penguin operators 3 - G , see the review [5] for precise definitions. 
For the considerations in this paper, one can safely neglect contributions of 
the penguin operators and consider only contributions of the leading Oi y2 
operators. The 0\ operator can be projected to a color-singlet state as 



1 = —0 2 + 20 2 



where 6 2 = (cT^(l - 75 )c)(sTV(l - 7s)&)- 

For H c = i] c , J/ip, Xco, Xc\i the decay rate can then be written as 



F(B -> KH C ) 



G 2 



F 



32tt 



v cb \ 2 \v:f 



i 

m B 



1 - 



rn 



m, 



\M Hc 



and 



(5) 



(6) 



M Hc =M n H onfact + M J H a 



.fact 



cm + 



CM 



(H c K\0 2 \B) + 2C 1 (n)(H c K\6 2 \B) 



(7) 



The first part of (7) can be calculated by factorizing the matrix element of 
the 2 operator as 



3 



(H K\O 2 \B) = {H c \(cc) V - A \0){K\(sb) V - A \B) , 



(8) 



and using the corresponding expressions for the (H c \(cc)v-a\0) from Table 
2. Note that due to the reason stated in introduction, Eq.(l), there is no 
factorizable contribution to the B — > K\ c q decay. The B — ^ K matrix element 
from (8) is defined by the decomposition 



(K(q)\s^b\B(p + q)): 
(2q + p) fM F+ K (p 2 ) + 



(9) 



while the estimation of the form factors F£ K and F BK in the LCSR approach 
[10] gives the following values needed in our calculation: 



f$K(™%+) = 0.60 , F+ K (m 2 x J = 0.74 , F° BK (m 2 J = 0.42 , 



(10) 



with the theoretical uncertainty of 15%. 

For a particular charmonium considered, the expression (7) can be then brought 
into the following form: 



M 



j/i> 



2e • Pk fj/^rnj/^F^ 



BK \ m J/i> 



M Xcl = 

■ Pk f Xcl m Xcl F£ K (m 2 Xci ) 



M Vc =tm 2 B f Vc F° K (m 2 c ) 
M Xc0 = 2C 1 (ti)m 2 B f Xc0 F x M- 



FMm 2 x J\ 



'ID 



Fh c =(ti c ,j/iI), Xc o, Xc1 ) represents the nonfactorizable part directly proportional to 
the contribution of the 2 operator in (7). This contribution vanishes under 
the factorization assumption. Below we present the calculation of the nonfac- 
torizable soft contributions in the LCSR approach. 
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2 Nonfactorizable contributions in B — > fT-f/c (iJ c = J/ip, XcJ(j=o,i)) 
decays from LCSR 

The starting point of the calculation using the LCSR method is the correlation 
function, defined as 

T = i 2 J dx 4 J dy 4 e-^^^- fe )^(X(p^)|j ifc (y)(5 2 (0)j B (x)|0), (12) 

where p 2 K = m 2 K = 0, k 2 = and p 2 B = m 2 B ,p 2 Hc = m 2 Hc . The interpolating 
current of a B~ meson is given as jb = ifnjy^^u, whereas the choice of the 
interpolating charmonium current has to be done according to the definite 
J, P, and C quantum numbers of a particular meson H c . The considered 
charmonium currents, together with some properties of charmonium needed in 
the sum rule analysis are summarized in Table 2. For the consistency, we also 
include discussion on the B — > KJ/ip decay which was already extensively 
presented in [4] and will only quote here the numerical result. For the rest 
of the B decays into charmonium, the calculation closely follows the LCSR 
approach developed in [3] and [4] and we refer to these references for all details. 
Including the twist-3 and twist-4 nonfactorizable contributions, we can write 
first the nonfactorizable contribution to the B — > K\ c i decay analogously to 
the result for B -> KJ/ip, Eq.(65) in [4]: 




where s c and Sq are the effective threshold parameters of the perturbative 
continuum in the H c and B channel, respectively, and Uq = (ml—m 2 Hc )/ (s B — 
m^j), here specified for H c = Xci charmonium. M, the Borel parameter in 
the B channel, and the parameter n = 1,2,... in the charmonium channel 
have to be chosen in such a way that a reliable perturbative calculation is 
possible, but on the other hand that excited and continuum states in a given 
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channel are suppressed. The function X appearing above is X = x(s, y, m B ) = 
(s — Am 2 / (1 — y))/(s — m B ) and the expansion up to 0(X 2 ) is performed. The 
twist-3, fox, and the twist-4, <p±, three-particle kaon distribution amplitudes 
are defined as usual [11]. The scale at which F H Js are calculated is fi b ~ 
m b /2 ~ 2.4 GeV. 

As for the B decays into the scalar and pseudoscalar charmonium the calcu- 
lation yields 



1 



1 du 



2rn 2 B f B J u b u 



^ e (m|-K-m2 c (l-„))/„)/M2 



1-- 



dy 



m b m c 



{i-y)Vv m vc 



i i ■ 

mi — m^ c 



— 03x(l -u,u-v,v) 

V \ u 



ru dy _ 

+3m b f K — ± (1 - u, u - v, v) 
Jo v 



(14) 



and 



ds 



>L + Ql) 



2\n+l 



L_ f 1 ^ e (m% 
r/ B Jug U 

2 ™2 



2%/b - -o 



-«))/«)/M 2 «" < 

7o (l- 



m b m c 



y)Vy m xc 

X\ 



f u dv , . (ml - m 2 n \ / , X 

Jo V \ U J \ V 



+3m b f K 



— 0_l(1 - u, u - v, v) (y + (1 - y)— ) 



(15) 



3 Numerical predictions and discussions 



Let us first specify the numerical values of the needed parameters. For param- 
eters in the B channel we use m B = 5.28 GeV and the values taken from [12]: 
f B = 180 ±30 MeV, m b = 4.7 ±0.1 GeV, and = 35 ±2 GeV 2 . For the char- 
monium states we use the parameters from Table 2 and m c = 1.25 ±0.10. The 
K meson decay constant is taken as fx = 0.16 GeV. For parameters which 
enter the coefficients of the twist-3 and twist-4 kaon wave functions we sup- 
pose that ,/W ~ f 3K and 5 2 K ~ 5%, and take f 3K = 0.0026 GeV, 5 2 (/i b ) = 0.17 
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Table 3 

Theoretical results for the B — > -fT-f/c decays calculated in this paper. FH c {^b) is the 
nonfactorizable contribution, M™™ * act /M.^* , Eq.(7), is the ratio of the nonfactor- 
izable and factorizable amplitudes for a particular mode B — > KH C , whereas B is 
the branching ratio, all calculated at ji = [if,- Large scale-dependent uncertainties 
pertinent to the factorizable amplitude are not included. 



decay mode 




M ™fact/ M fact 


s(io- 4 ) 


B~ -► K-r] c 


0.0015-0.0019 


0.08-0.10 


2.0 ±0.1 


B~ -► K~Jip 


0.011-0.018 [4] 


0.40-0.70 [4] 


3.3 ±0.6 


B~ ^"Xco 


-(0.0007-0.0008) 




0.0017 ± 0.0002 




0.044 - 0.052 


1.30-1.50 


5.1 ±0.5 



GeV [11]. The stability region for the Borel parameter is found in the interval 
M 2 — 10 ± 2 GeV 2 , known also from other LCSR calculations of B meson 
properties. Concerning the sum rules in the charmonium channels, the cal- 
culation is rather stable on the change of n in the interval n = 4 — 7. Qq 
is parameterized by Ql = 'im 2 ^, where in order that the lowest resonances 
dominate £ takes values from 0.5 to 1 for the B decays into s-wave charmonia, 
while £ is between 1 and 2.5 in the decays into p-wave charmonia [7]. 

The results of our calculation are summarized in Table 3. Comparing the 
numbers from Table 1 and Table 3, it is important to note that in general 
the calculated branching ratios are still too low to accommodate the data, 
except maybe for the B — > K\ c i decay. The nonfactorizable correction to 
B — > KJ/ip and particularly to B — > K\c\ is large, for the B — > K\ c i decay 
this correction is even larger than the factorizable contribution, Table 3. On 
the other hand, in B decays into (pseudo) scalar charmonia the nonfactoriz- 
able contributions are small. The reason is the cancellation of the twist-3 and 
twist-4 contributions in these decays. However, even without this cancellation, 
the nonfactorizable effects produced by the exchange of a soft gluon between 
a kaon and Xco as calculated here, could not be able to account for such a 
large branching ratio of B — * Kxco as measured experimentally. These would 
demand contributions which have to be at least an order of magnitude larger 
than those typically occurring in a B decay into charmonium state as esti- 
mated by the LCSR method. Therefore it is very unlikely that the mechanism 
of the nonfactorizable soft gluon exchange is the reason for a relatively large 
branching ratio of the B — > Kxco decay. A possible speculative explanation 
for large discrepancies between the theoretical results for the B decays into 
(pseudo) scalar charmonia and the experimental data one could find in the 
nonperturbative effects of the instanton type. They could appear due to the 
light quark admixtures in the i] c and Xco mesons, yet it would be hard to 
account for such contributions reliably. 
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The B — > Kxco decay was also analyzed within QCD factorization method [13] 
and it was observed [14] that there is a problem of logarithmic divergences of 
the decay amplitude already at the leading-twist order. In another approach, 
by studying the mechanism of the rescattering of charmed intermediate states 
in B — > Kxco, the authors of Ref.[15] show that the rescattering effects could 
provide the large part of the B — > Kxco amplitude. 

To summarize, we have studied the soft nonfactorizable contributions to B — > 
KH C (H C = rj c , J ftp, Xco, Xci) decays by using the LCSR approach. Inspite of 
the expected large contributions which could explain large discrepancy be- 
tween the factorizable predictions and the experimental data, we were not 
able to confirm these expectations, except for the B — > Kxci decay, for which 
large nonfactorizable corrections are found. The other predicted B decays 
into charmonium receive nonfactorizable soft contributions too small to ac- 
commodate the data. Unfortunately, this is particularly true for the puzzling 
B — > Kxco decay that factorized amplitude vanishes and the LCSR mech- 
anism considered in this paper cannot explain its relatively large branching 
ratio. 

Note added: 

While this work has been prepared, Ref. [16] appeared, where the authors 
discuss nonfactorizable soft contributions in the B — > Kr\ c and B — > Kxco 
decays within the LCSR approach. There are some differences observed in 
the approach as well as in the results between [16] and our paper. As for the 
B — > Ki] c decay, the authors of [16] choose the pseudovector current for i] c . 
In that case, the result can be easily derived from Eq.(13) according to the 
approach taken from [4]. The differences show in the twist-4 part. Apart from 
this, by taking the pseudovector current to describe r\ c one has also to include 
the mixing with the Xci(l ++ ) state explicitly [17] which was not considered in 
[16]. More importantly, for the problematic B — > Kxco decay, apart from the 
superfluous factor x in Eq.(49) of [16], we agree analytically in the twist-3 part. 
In the second version of [16] , the authors have corrected numerical errors and 
have included the twist-4 contribution to their result for the nonfactorizable 
B — > Kxco amplitude. Although we now agree that the soft nonfactorizable 
contribution in the B — > Kxco decay is small and cannot accommodate the 
experimental data, we still disagree in the twist-4 part of the nonfactorizable 
contributions, which renders the numerical results from this paper somewhat 
different from those obtained in [16]. 
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